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ABSTRACT 

We present a magnitude and proper motion limited catalogue of ~ 10, 000 white 
dwarf candidates, obtained from the SuperCOSMOS Sky Survey by means of reduced 
proper motion selection. This catalogue extends to magnitudes R ~ 19.75 and proper 
motions as low as ~ 0'.'05 yr""'^, and covers nearly three quarters of the sky. Pho- 
tometric parallaxes provide distance estimates accurate to ^ 50%. This catalogue is 
used to measure the luminosity functions for disk and spheroid white dwarfs, using 
strict velocity cuts to isolate subsamples belonging to each population. Disk luminos- 
ity functions measured in this manner are really a conglomerate of thin and thick disk 
objects, due to the significant velocity overlap between these populations. We intro- 
duce a new statistical approach to the stellar luminosity function for nearby objects 
that succesfuUy untangles the contributions from the different kinematic populations, 
without the need for stringent velocity cuts. This improves the statistical power by 
allowing all stars to contribute to the luminosity function, even at tangential veloci- 
ties where the populations are indistinguishable. This method is particularly suited to 
white dwarfs, for which population discrimination by chemical tagging is not possible. 
We use this technique to obtain the first measurement of the thick disk white dwarf lu- 
minosity function, while also improving constraint on both the thin disk and spheroid 
luminosity functions. We find that the thin disk, thick disk and spheroid populations 
contribute to the local white dwarf density in roughly 79%/16%/5% proportions. 

Key words: White dwarfs; stars: luminosity function, mass function; surveys; Solar 
neighbourhood. 



1 INTRODUCTION 

Nearly all stars in the Galaxy will end their lives as faint 
white dwarfs (WDs), once nuclear burning has ceased and 
the outer photospheric layers have been ejected. The further 
evolution of these objects is a gradual cooling process, where 
the remaining thermal energy of the exposed core is slowly 
radiated away through the thin outer layers. WDs are highly 
numerous in the Solar neighbourhood. They are the second 
most common type of star behind low mass main sequence 
stars, but due to their extreme faintness are found in much 
smaller numbers in Galactic star surveys. As a result, the 
evolution of the coolest WD s is rather poorly constrained 
observationally (iHansen fc L icbcrt 2003; Gates gt al. 2 0041 : 
IVidrih et al.ll2007l : lRowell et al.ll2008i : iKilic et al.ll2010l ). and 
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though much work has been done on the theoretical front in 
recent years (|Fontaine et al.l 1200 ll : [Bergeron et al.l 1200 ll ) it 
is still unclear how to interpret t he small amount of obser- 
vational evidence available (e.g. iBergeron fc Leggettll200^ : 
iBedin et al.l |2008| ) . Studies of the luminosity function for 
white dwarfs (WDLF) typically focus on the faint end, where 
the finite age of the Galaxy predicts a sudden downturn in 
the local density of WDs when the cooling times converge 
on the age of the popula tion as a who l e. Th e first evidence 
for this was published b y lLieberT et al.l(|l988l ). and was used 
bv lWinget etal] l| 19871 ) to obtain an estimate of the age of 
the universe essentially independent of main sequence stel- 
lar ages. Although many more studies have been carried out 
in the i ntervening d ecades, focu ssing on both the Galac- 
tic di sk l|Oswah et a l. 1996; Lcgg ett et al][l99 8: Kn ox et all 
[1999') and clusters (Hansen et al.ll2002l : [Bedin et al.il2010l ') it 
is the advent of digitized wide angle sky surveys in the last 
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decade that has revolutionised the study of these objects in 
the field, with large footprint areas greatly increasing the 
survey volume for faint objects. 

The Sloan Digital Sky Survey (SDSS) has provided the 
most interesting results in recent years. Several indepen- 
dent surveys have increased the catalogue of known WDs 
by more than an order of magnitude, providing greater sta- 
tistical power at both the hot end of the LF, where the 
WDs are intrinsically rare, and the cool end where they ex- 
ist in large number bu t are extremely faint. For example, 
iKrzesinski et al.l l|2009h use colour selection to identify a few 
hundred hot WDs from the SDSS spectral catalogue. Their 
luminosity function extends to Mi,oi = 7.0, where the sample 
starts to become incomplete due to the selection probability 
for spe ctroscopic follow up dr opping off at lower temper- 
atures. lEisenstein" et al.l (|2006l ) identified more than 9,000 
spectroscopically confir med WDs in SDSS DR4, and their 
catalogue was used by iDe Gennaro et all (j2008l ) to probe 
the WDLF for different WD mass ranges. As cool WDs are 
generally selected for SDSS spectroscopic follow-up by algo- 
rithms designed to target QSOs, the complete ness of such a 
cat alogue is a strong f unction of colour. Kilic et al.l (|2006l ) 
and lHarris et al.l (j2006l ) (hereafter H06) used reduced proper 
motion to identify ~ 6000 WDs, with proper motions de- 
termined by pairing sources with the USNO-B catalogue. 
The reduced proper motion approach offers a superior sur- 
vey method for a couple of reasons. First, it is insensitive 
to colour, and cool WDs can be identified among the highly 
numerous disk stars from which they are indistiguishable in 
colour space. Second, it provides a statistically rigorous se- 
lection probability, avoiding the need to calculate detailed 
completeness functions to prevent survey bias. The WDLF 
of HOG is currently the deepest and best constrained around 
the faint turnover in the disk LF, though the pairing cri- 
terion limits the depth of their catalogue to the magnitude 
limits of the USNO-B survey. 

There has been considerable interest in Galactic 
spheroid WDs in the last decade, though very little is known 
about this elusive population of stars. Broad predictions can 
be made on the bas is of Galactic formation and stellar evo- 
lution models, (e.g. lHansenlll999l : lReidll2005h but the contri- 
bution made by these objects to the local density and Galac- 
tic mass budget is very poorly constrained. This is due to a 
combination of several factors. Firstly, these objects are both 
faint and intrinsically very rare, and in any given survey will 
be found in much smaller numbers than both disk WDs and 
other types of star. They are far too faint to be seen at large 
Galactic plane distances in regions of pure spheroid stars, 
and any objects must be observed as they pass through the 
nearby Galactic disk. They also cannot be identified by low 
metallicities, because for all WDs the surface gravity is high 
enough to pull metals below the photosphere. The only reli- 
able way to identify them is by their large space velocities, 
as the high velocity tail of the distributi on is cleanly sepa- 
rated from both the thin and thick disk (IVidrih et al ] |2007l : 
iKilic et aLllioiol ). This technique was used by HOG to mea- 
sure the first luminosity function for spheroid WDs of any 
significant size. 

The thin and thick disk WD populations overlap too 
much in velocity for single objects to be reliably attributed 
to either population. For this reason, there has never been 
a direct measurement of the thick disk WDLF, though all 



previous measurements of the WDLF for thin disk stars will 
have contained a significant number of these objects. In this 
paper, we show that the thin and thick disk WDLFs can 
indeed be measured separately from a mixed catalogue, by 
avoiding the need to classify each survey object as one type 
or the other. This technique has more statistical power than 
the traditional approach, because all objects are allowed to 
contribute to the luminosity function, even at tangential ve- 
locities where the populations are indistinguishable. 

This paper is organised as follows: in Section [2] we de- 
scribe the source data and detail efforts to construct a com- 
plete catalogue of stars. In Section [3] we use proper motion 
survey techniques to obtain a subsample of WD candidates 
from our catalogue, and in Section |4] we review the cata- 
logue with emphasis on the completeness and likely mixture 
of kinematic populations present. In Section [S] we describe 
the traditional — approach to measuring the luminosity 
function from a proper motion catalogue, and in Section [S] 
we use this technique to measure luminosity functions for 
various velocity subsamples drawn from our WD catalogue. 
In Section [7| we present our statistical approach to separat- 
ing the kinematic populations, and measure WDLFs for the 
thin disk, thick disk and spheroid. In Sections [8] and |9] we 
compare our results to those of other studies and draw our 
conclusions. 



2 THE DATA 

The SuperCOSMOS Sky Survey (SSS) was compiled by dig- 
itizing several generations of photographic Schmidt plate 
surveys. The source material includes the second genera- 
tion POSSU-B,R,I and SERC-J/ EJ, ER/AAO-R, I sur- 
veys, observed by the Palomar Oschin Schmidt Telescope in 
the northern hemisphere and the United Kingdom Schmidt 
Telescope in the south. These were performed in the photo- 
graphic bj r-59F ijv bands. The subscripts refer to particular 
combinations of filter and photographic emulsion used by 
the surveys; we ignore the small differences in response be- 
tween the northern and southern hemisphere. 

The first generation POSSI-i5 and ESO-R surveys pro- 
vide an early epoch for astrometric constraint. The photom- 
etry is mostly redundant due to a similar response to r^gp, 
but we will nevertheless refer to magnitudes in these surveys 
as rio3aE and resF when distinction has to be made. Note 
that rio3aE and resF are identical in response, but we do 
expect differences in the survey characteristics due to e.g. a 
different copying history prior to digitization. These surveys 
cross over at S ~ —18°, thus the sky is divided into three 
regions of common photographic material. This has impli- 
cations for the charcteristics of the digitized survey. Over 
1700 fields are required to cover the whole sky, with each 
field being observed in each of the four photographic bands 
to a depth of r ~ 21. The four observations in each field 
are spread over a wide time baseline in the latter half of the 
twentieth century, with up to fifty years between epochs in 
the extreme. 

Schmidt plates were digitized by the SuperCOSMOS 
scanning machine, scanning at 0.67 arcsecond pixel size 
and 15 bit digitization. SuperCOSMOS was operated at the 
Royal Observatory of Edinburgh - see Hambly et al. (2001a) 
for an introduction to SuperCOSMOS and brief historical 
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census of scanning programs. In some cases, glass or film 
copies of survey plates were scanned rather than originals. 
This includes the entire northern hemisphere. The use of 
photographic copies was not thought to significantly degrade 
the quality of the scanned data, however we regard each of 
the surveys as having independent characteristics in case any 
differences arise. Image Analysis Mode software was used to 
process the raw pixel data into parameterised object cata- 
logues, which were then merged across the four epochs to 
measure proper motion and colours. 

The SSS data are housed in a relational database 
at the Wide Field Astronomy Unit, Royal Observa- 
tory of Edinburgh. The original SSS is now contained 
within the SuperCOSMOS Science Archive, with ac- 
cess provided via a web-based SQL query form at 
http://surveys.roe.ac.uk/ssa/. Parametrised object in- 
formation is dispensed from two large tables; the detection 
table contains roughly 6.4 billion individual object detec- 
tions, which are merged into around 1.9 billion multi-colour, 
multi-epoch observations in the source table. A comprehen- 
sive description, user-guide and technical notes for the SSS 
were released in a series of three papers - see Hambly et al. 



2.1 Proper motion limits 

The input data for our proper motion survey is drawn 
from two sources. The online SSA provides low proper mo- 
tion coverage, from a rough lower 5a detection limit of 
CK'OS yr~^ to a fixed upper limit of O'.'lS yr~^. The com- 
pleteness of the SSA at this upper limit varies widely be- 
tween fields, due to the differing epoch spread and fixed 
search radius for object pairing. Early tests revealed that 
fields start to become incomplete above around C/.'OS yr~^, 
as indicated by a significant departure from the linear rela- 
tion log(EiV(> /x)) oc log(^). Our low proper motion data 
is thus truncated at an upper limit of O'.'OS yr~^, with lower 
limits set according to the procedure outlined in the follow- 
ing section. 

Faster moving objects are recovered usi ng a more so- 
phisti cated search algorithm described in iHamblv et al.l 
l|2004h . Object pairing between epochs works on a field-by- 
field basis, and starts with the complete set of parameterized 
object detections on each of the four plates. Any sources that 
have been succesfuUy merged and included in the default 
SSS catalogue at motions lower than the O'.'OS yr~^ limit are 
thrown away. The remaining objects are then paired in all 
possible combinations within a search radius set by the cho- 
sen upper proper motion limit of lO'.'O yr~^ and the epoch 
difference between plates. The primary pairing is between 
the two r epochs, which are subject to a magnitude limit of 
'"59F < 19.75 to reduce noise, but any detections in bj and 
In are folded into the analysis to improve the astrometric 
solution. To produce an all sky catalogue, we take the entire 
object catalogue and purge multiple observations appearing 
in plate overlap regions, keeping whichever appears closest 
to its respective field centre in order to match the seaming of 
the low proper motion catalogue. This is far more rigorous 
than the object pairing in the default SSS, and is designed to 
maximise completeness of high proper motion objects at the 
expense of introducing large numbers of spurious detections. 



Appropriate selection of astrometric and image statistics is 
necessary to reduce the contamination to a tolerable level. 



2.1.1 Low proper motion completeness limits 

We impose a lower proper motion limit on a field-by-field 
basis such that all objects have at least a 5a^ proper motion 
detection, defined by 
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This excludes non-moving objects from our catalogue, and 
limits scatter in the reduced proper motion diagram that 
is used to select white dwarf candidates. However, we can- 
not simply select all objects with ^ 5(t,j detections, because 
the proper motion errors show significant scatter at con- 
stant magnitude. The resulting survey volume limits at given 
tangential velocity would be unknown, as it would be im- 
possible to measure the distance at which the star dropped 
below the required Sct^ detection threshold. Therefore, we 
wish to find the maximum proper motion error as a func- 
tion of magnitude, and use this to fix the lower proper mo- 
tion limit. This guarantees that all objects that pass the 
limit also have ^ 5(7^ detections. The proper motion er- 
ror varies from field to field, due to differences in plate 
quality and time baseline. It is also a strong function of 
magnitude, and at constant magnitude shows a significant 



spread. Figure 1(a) shows the distribution of cr^ with bj for 



a representative field. We use bj in this analysis because it 
is the highest quality photometry available. Objects in the 
upper locus have been missed at one of the four epochs, 
usually rasF , and have inferior astrometric fits. For this rea- 
son we restrict our low proper motion sample to objects 
with four plate detections. To fix the proper motion limit, 
we fit functions to the upper boundaries of the 6j,o"^ locus. 
The procedure for each field is as follows. We start at the 
bright end of the distribution, and bin objects on magni- 
tude interactively so that all bins contain one hundred ob- 
jects. The mean bj and maximum in each bin is located, 
after rejecting the top 5% of ct^ as outliers. The full set of 
(bj) ,0-™"^ points defines the rough upper boundary of the lo- 
cus, but shows considerable noise on small scales. The raw fit 
is then processed through one stage of smoothing to obtain 
the desired function. We use a Savitzky-Golay technique to 
smooth the data, general ized from their original prescription 
(|Savitzkv fc Golavlll964l ) to allow for non-uniform points on 
the abscissa. This removes small-scale noise while preserv- 
ing any low-frequency features, such as the mid-magnitude 
turning point seen in most fields. Figure [T(b)| shows the same 
field as (a), but restricted to four plate detections, and with 
the fitted functions shown. The final, smoothed function ob- 
tained for each field is denoted a'^°'^{bj), and is used to set 
the lower proper motion limit according to 

M™n(6,/) = 5(a™"(6j) -f 0'.'002 yr'^). 

The set of points defining o-™"^(&j) is tabulated for each 
field, and interpolated using cubic splines to obtain the 
proper motion limit at arbitrary apparent magnitude. Thus 
the lower proper motion limit is a piecewise function of ap- 
parent magnitude, which complicates the measurement of 
the survey volume, but maximises the sample size by fully 
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Figure 1. Field 772S proper motion errors and fitted o-^{bj) 
function, (a) includes objects not detected at rgsFi which form 
the locus of points at higher o-^. (b) is restricted to objects with 
full four plate detections, and shows the function fitted to the 
upper boundary of the (7^,6 j locus. 



2.1.2 High proper motion 

The precise location of the upper proper motion limit is 
not as important as the lower proper motion limit from a 
catalogue completeness point of view, because of the low 
probability of finding objects at these velocities and the fact 
that the survey volume is not as sensitive to errors at high 
proper motions {ay oc n~^a^). However, we expect to find a 
greater fraction of elusive spheroid WDs at high proper mo- 
tions, and extending the range as wide as possible while lim- 
iting contamination from spurious detections allows greater 
constraint on the luminosity function. 

With far fewer high proper motion objects present, 
there is insufficient statistical power to set an upper limit on 
a field- by-field basis and a global approach is required. We 
investigated the completeness of proper motions by looking 
at the cumulative number counts for all objects, after ap- 
plying the selection cuts outlined in the remainder of this 
section to eliminate noise. Figure [3] shows the cumulative 
number counts as a function of proper motion for all sur- 
vey objects, extending into the low proper motion range to 
demonstrate the completeness across the interface as indi- 
cated by the first vertical bar. We fix the upper limit at the 
point where the number counts depart from a straight line, 
which was judged to he at fj, = I'.'O yr~^. This is marked by 
the second vertical bar. 



exploiting the variable proper motion accuracy. The mean 
lower proper motion limit across the survey is shown in Fig- 
ure [2] This is divided into the three sky regions where the 
source photographic data are the same. Note that the lower 
quality northern data attains a similar proper motion accu- 
racy as the southern region that shares the same first epoch r 
material, that of the POSSI-E survey {S >-18°). The ESO-R 
survey has a much later average epoch, with the result that 
proper motion measurements are considerably more uncer- 
tain. 



2.1.3 Astrometric residuals 

We place restrictions on the astrometric xl in order to ex- 
clude objects with a poor fit to their proper motion solution. 
Appropriate cuts are obtained iteratively by checking for 
contamination in the reduced proper motion diagram and 
the distribution of for the resulting catalogue. A cut 

of xi ^ 2.0 is applied to the high proper motion data across 
the whole sky, and xt ^ 1-5 in the southern hemisphere for 
the low proper motion data. The low proper motion data 
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Figure 3. Completeness limit at high proper motion. 



each photographic survey, b and f! are the Galactic latitude and 
solid angle subtended by each field. 



is adversely affected in the north by the lower quality plate 
material, and a restrictive cut of xt ^ 0.8 is necessary. 

The stars ejected from the survey by these cuts can 
be accounted for later by including an appropriate correc- 
tion factor in the LF measurement. In both the high and 
low proper motion data, the requirement of a detection at 
all four epochs results in a four degree of freedom astro- 
metric fit {u = 4). However, the distribution of xt for ob- 
jects in both the SSA and our high proper motion data is 
significantly tighter than expected. We have analyzed the 
astrometric residuals for objects in a selected SSS field in 
conjunction with synthetic data, and established that the 
tangent-plane position errors used in the astrometric fits 
are overestimated by a factor of around \/2, leading to a xt 
roughly half its true value. This makes the correction factor 
hard to determine accurately. Assuming all xi values are 
halved, & xl ^ 2.0/1.5/0.8 cut will eject around 0.5/2/20% 
of stars. 



2.2 Magnitude limits 

The SSS literature provides magnitude completeness esti- 
mates for a selection of bj and r^gp plates in the first SSS 
data release, the South Galactic Cap (SGC). These are mea- 
sured by comparing SSS star and galaxy counts with those 
obtained from deeper prime-focus and CCD data in overlap- 
ping regions, and estimate near 100% completeness within 
~ 1.5 mags of the plate detection limits. However, our all- 
sky survey uses northern Schmidt plate data not included in 
the SGC, those of the POSS-I and POSS-H surveys. Also, 
the r63F/i03aB and in plates are shallower and of lower sig- 
nal to noise than bj and rsgp, and, given the colours of the 
objects we are interested in, will likely determine the overall 
completeness limits of the SSS fields when considering ob- 
jects detected on all four plates. Therefore, it is necessary to 
estimate new completeness limits for all plates used in the 
SSS. Since we cannot obtain dee per imaging for all fields, 
we instead follow the method of iTinnev et al.l l|l993h : this 
involves simulating star and galaxy counts along the line 
of sight, and comparing these to observed counts derived 
from the corresponding plate material. We assume that the 
completeness characteristics of plates within the same pho- 
tographic survey are identical, due to uniform quality con- 
trol, emulsion grade and copying history prior to digitiza- 



tion. This allows us to restrict our analysis to a represen- 
tative sample of plates from each of the eight photographic 
surveys used in the SSS. Although the plate detection limit 
varies within each survey, we assume that the plates have a 
common completeness function, which we define as the ratio 
of detected objects to real objects as a function of magni- 
tude relative to the plate detection limit. We analyzed five 
plates from each survey, drawn from five fields in each of the 
celestial hemispheres. A summary of the fields selected for 
analysis is given in Table [T] 



2.2.1 Synthetic star and galaxy counts 

2.2.1.1 Stars Differential star counts along the line of 
sight to eac h field are obtaine d using the Besangon Galaxy 
model (see I Robin et al.l '2003) . This employs a population 
synthesis approach to produce a self-consistent model of the 
Galactic stellar populations, which can be 'observed' to ob- 
tain theoretical data sets useful for testing various Galactic 
structure and formation scenarios. We use the coordinates, 
solid angles and passbands of our selected fields as inputs, 
and select an output number count range that goes several 
magnitudes fainter than the plate limits. The SuperCOS- 
MOS filter system is not included in the Besangon model; 
instead we use the CFHT Megacam bands g and r to ap- 
proximate bj and r^gp^ggp^iQgaE , and Johnson-Cousins / to 
approximate ijv. 

2.2.1.2 Galaxies Within a few magnitudes of the plate 
limits, galaxies appear as unresolved, point-like objects and 
have image parameters that overlap with stars. We therefore 
have to include galaxies in our synthetic number counts. 
Galaxy counts to faint magnitudes have been determined 
in many independent studies. We use counts provided by 
the Durham C osmology Group t hat combine their ow n re- 
suhs (see e.g. Ijones et al1ll99ll : [Metcalfe et al.lll99lh with 
those of many other authors. These are available onlin^B 
along with transformations to photographic bands. They 
are provided in terms of log-number counts per square de- 
gree per half-magnitude; we fit straight lines to obtain func- 
tional forms for the galaxy counts in each band, and trans- 
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form these to O.IM for comparison to our observed counts. 
Note that the ~ 25 square degree field of view of each 
Schmidt plate smooths out any anisotropics in the faint 
galaxy counts. The fitted functions are given in Equations [1] 
(r = r59j7/g3j7/io3a£;). The units on the number counts A'' are 
deg~^ 0.1M~^. We multiply these functions by the solid an- 
gle of each field, then add them to the star counts to obtain 
the total synthetic counts for each plate. 

\og{Nb,) = 0.471 bj - 7.890 (16 < h., < 26) (1) 
log(iVO = 0.379 r - 5.351 (17 < r < 25) 

log(iVi„) = 0.606 iiv - 9.132 (12 < ijv < 17.75) 

= 0.346 iiv - 4.397 (17.75 <iN < 25) 



2.2.2 Observed star and galaxy counts 

For each field selected for analysis, we obtain object counts 
to the detection limit on all four plates, binned at 0.1 mag- 
nitude intervals server-side using an efficient SQL query on 
the SSA interface. We use a wide range of image statistics to 
include partially resolved galaxies, though close to the plate 
limits everything is pointlike. 
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Table 2. Offsets defining the magnitude limit for each field in the 
eight constituent photographic surveys. Column three gives the 
average survey magnitude limit on applying these offsets globally. 
The number in brackets is the corresponding average when the 
''59F < 19.75 constraint is applied to the second epoch r data, as 
explained in Section 12.2.41 



large offset of 1.9m. This is a significant restriction on the 
survey volume, considering that the POSSI-i5 survey covers 
the whole sky down to 5 ~ —18°. 



2.2.3 Completeness functions and famt magnitude limits 

In Figure |3] we show observed and modelled differential ob- 
ject counts for for each plate in field 270 in the north. The 
ratio of these quantities relative to the detection limit gives 
the completeness function for the plate; these are inlaid for 
comparison. Total model counts are normalised to the ob- 
served number at two magnitudes above the detection limit, 
where the plates are assumed 100% complete. To estimate 
the global completeness function for each survey, we repeat 
the analysis for the five representative fields listed in Table[T] 
and take an unweighted average of the individual complete- 
ness functions. The global completeness functions obtained 
for the northern hemisphere surveys are shown in Figure [5] 
for reference. We set the magnitude limits in terms of an off- 
set from the detection limit, defined as the magnitude of the 
faintest detected object. This ensures that the noise at faint 
magnitudes is avoided. The global completeness functions 
are used to judge an appropriate value. The offsets selected 
are listed in Table (2] along with the mean magnitude limit 
obtained on applying the offset to all survey plates. 

The superior quality of SERC-J and SERC-ii is evi- 
dent, due to the use of original glass survey plates in the Su- 
perCOSMOS scanning program. These plates show ~ 100% 
completeness to within a few tenths of a magnitude of the 
detection limit. All other surveys were copied photograph- 
ically at least once before digitizing, which has resulted in 
noise creeping in within a magnitude or so of the detection 
limit. The POSSll-J survey was copied twice before scan- 
ning, and the noise is noticeably worse on these plates. The 
POSSF-E plates often show extremely large numbers of spu- 
rious detections within up to two magnitudes of the detec- 
tion limit. A double peak is often seen in the distribution; 
this is due to the mosaicking of two or more Palomar fields 
of different depths onto one ESO-SERC field for inclusion in 
the SSS catalogue. This means that to ensure completeness, 
the POSSI--E magnitudes have to be limited to the depth 
of the shallowest plate in mosaicked fields, leading to the 



2.2.4 Astrometric noise at faint magnitudes 

Early catalogues drawn using the magnitude limits de- 
scribed in the preceeding Section showed large amounts of 
contamination by faint objects with erroneous proper mo- 
tions. This is most apparent in the reduced proper motion 
diagram that is used to select white dwarf candidates (see 
Section [Sjl, where large numbers of faint disk main sequence 
stars with erroneous large proper motions scatter into the 
region populated by cool white dwarfs. Many attempts were 
made to restrict these objects using astrometric and im- 
age statistics, but in the end a survey-wide cut on the sec- 
ond epoch r69_F magnitude was required. We fix this cut at 
r 19.75. r59F is used in the primary multi-epoch object 
pairing, so is a good proxy for the robustness of the proper 
motion measurement. The numbers in parenthesis in Table 
[2] give the mean magnitude limit on combining this cut with 
the completeness limits. 

This involves modelling the effective survey volume over 
discrete tangential velocity ranges, and turns out to be a 
good way to separate the disk and spheroid white dwarf 
populations, which are distinguishable only on the basis of 
kinematic study. 

2.2.5 Bright magnitude limits 

The selection of the bright limits on apparent magnitude is 
not as crucial, due to the low probability of white dwarfs 
being found at these magnitudes. We fix the bright limits 
at 12 for all bands and across all fields. This is within the 
complete range of the Schmidt plates. 

2.3 Photometric accuracy and transmission 
functions 

While photography cannot compete with CCD astronomy 
in terms of photometric accuracy, it has been noted in the 
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Example completeness analysis for two fields in the SSS, covering all eight of the constituent photographic 



literature (e.g. ISalim et al.ll2003 ) that insofar as digitized 
photographic surveys go, SSS photometry is of the highest 
quality with uncertainties as low as am ~ 0.07. Rough esti- 
mates of the external error on the photographic colours at 
bj ~ 16.5 and bj ~ 20 are given in Hambly et al. 2001b 
These are used when fitting photometric models to stellar 
colours. We fit a straight line for the error at intermediate 
and fainter magnitudes, and use the uncertainty at bj ~ 16.5 
for all brighter magnitudes. The relation we obtain is 

,^ , fo.Q7 if fej < 16.5 

ffb-r b-i(bj) = S 

[O.Q26fej - 0.35 otherwise. 

Due to the particular way in which the photometric scale 
is calibrated, colours are more accurate than single magni- 
tudes. Uncertainty in single passbands is necessary however 
for deriving minimum-variance estimates of the photomet- 
ric distance, by comparing model and observed magnitudes. 
These are taken from Table 12 in Hambly et al. 2001b. Filter 
transmission fu nctions have b een obtained from from lEvand 
l|l989l ) (6j) and|Besseil (llOSel ) (resF/iosas, ^ssf and in)- 



2.4 Image quality criteria 

Every parameterized object detection in the SSA is accom- 
panied by a set of image statistics, some of which we restrict 
in order that stars included in our survey have high quality, 
stellar images. 



Blend number We reject objects that have been de- 
blended at any of the four epochs. The de-blending algo- 
rithm attempts to recover individual object parameters, but 
is known to be unreliable. 

Quality number The quality number is a 32-bit integer 
flag set so that increasingly significant bits indicate in- 
creasingly compromised situations encountered during im- 
age analysis. We restrict this parameter to values less than 
128, which indicates an image likely to be afi'ected by a step 
wedge or other label on the photographic plate. 

Profile classification statistic The profile classification 
statistic rj provides a magnitude-independent measure of the 
'stellarness' of each image, by quantifying the deviation of 
the radial profile slope from that of a mean stellar template. 
rj is given in terms of a unit Gaussian statistic, and we accept 
images with \ri\ ^ 4(t. 

Ellipticity Previous studies utilising digitized Schmidt 
plate data have placed cuts on the ellipticities e of images, in 
order to limit contamination from faint galaxies and noise. 
However, we have found that the ellipticities provided by 
the SSA are extremely noisy at intermediate to faint mag- 
nitudes, and that any intuitively sensible cut will result in a 
seriously incomplete sample of stars. For example, real stel- 
lar images in the bj band are often assigned e > 0.5 within 
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Figure 5. Completeness functions averaged over five plates in each photographic survey. 



two magnitudes of the plate detection limit. We have ignored 
this parameter altogether. 



2.5 Sky Coverage 

2. 5. 1 Individual field areas 

The pointings used for the ~ 1700 fields comprising the 
SSS use the ESO/SERC system of field centres, which is 
based on a 5° pitch angle plus some small adjustments to 
allow for locating guide stars. This results in ~ 0.5° overlap 
between neighbouring fields, given the ~ 6° x 6° field of 
view of Schmidt plates. Objects observed multiple times in 
overlap regions are assigned to the field whose centre they 
are closest to along a Great Circle, providing a 'seamless' 
catalogue. The coverage of each field is not given in the SSS 
database, and has been measured for this work. This was 
done numerically by dividing the sky into small elements of 
solid angle, accounting for any excluded regions close to the 
plane, and assigning each element to the field it lies closest 
to along a Great Circle. The average field of view for fields 
in the SSS is ~ 0.007 sr, with a significant spread as shown 
in Figure [6l 
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Figure 6. Frequency of non-overlapping SSS field areas. 



accounted for in order to accurately calculate the survey vol- 
ume. The drilling fraction is measured exclusively from the 
h,j plates, on which these images have the greatest areal ex- 
tent, and ranges from a few percent at high Galactic latitude 
to more than 50% in the plane. 



2.5.2 Bright star haloes and diffraction spikes 

The fraction of each survey plate lost to internal reflections 
and other spurious images can be significant, and must be 



2.5.3 Rejection of crowded and dusty fields 

The SSS nominally covers the entire sky. However, in dusty 
regions at low Galactic latitude the photometric calibration 
can be adversely affected by differential extinction across 



White Dwarfs in the SuperCOSMOS Sky Survey 9 



Survey 


n{sv) 


Fraction of whole sky 


Population 


(U) 


{V) 


{W) 


0"!7 


av 


aw 


H < O'.'IS yr-i 


9.11 


0.72 


Thin disc 


-8.62 


-20.04 


-7.10 


32.4 


23.0 


18.1 


tM > O'.'IS yr-i 


8.80 


0.70 


Thick disc 


-11.0 


-42.0 


-12.0 


50.0 


56.0 


34.0 








Halo 


-26 


-199 


-12 


141.0 


106.0 


94.0 



Table 3. Total sky coverage of the low and high proper motion 
ranges. 



the field. Also, in crowded regions a large amount of de- 
blending is required, resulting in significant incompleteness 
as these stars are rejected by our survey procedure, and 
source merging across multiple epochs suffers from serious 
confusion. For these reasons, we avoid the Galactic plane by 
10° and the Galactic centre by 20°. We also reject 53 fields 
lying slightly outside this region that show a large amount 
of contamination, manifest in highly dense regions of stars 
bounded by field edges. Seven fields centred on the cores of 
the Magellanic clouds are also excluded. 



2.5.4 Rejection of fields with a small epoch spread 

In most fields the earliest epoch of observation is rggp/iog^E, 
with the remaining three distributed over a ~ 10 to 20 year 
period. The primary object pairing is between the two r 
epochs, and if the bj and ijv observations are taken very 
close to the r^gp observation, they provide very little as- 
trometric constraint. In such large number of first- 

epoch pairings are generated for every real high proper mo- 
tion object. As the search algorithm proceeds through every 
possible combination within the search radius, the single 
correct image pairing is swamped by noise. Tests show that 
at proper motions lower than around O'.'IS yr~^ this is not 
a problem, as the search radius for image pairing is small. 
However, at larger proper motions the effect can be catas- 
trophic. We therefore reject objects travelling faster than 
O'.'IS yr~^ in any field that has bj, rsgp and iat taken within 
1.5 years. 5S fields fall into this category, and are excluded. 
This results in slightly different sky coverage depending on 
the proper motion range. 



2.5.5 Total sky coverage 

The total survey footprint at low and high proper motions 
is given in Table [S] This excludes the Galactic plane and 
centre regions, accounts for the bright star drilling fraction, 
and all fields rejected from each proper motion range. 



3 SURVEY SELECTION CRITERIA 

3.1 Reduced proper motion selection 

The proper motions of nearby stars correlate with distance, 
in the sense that closer objects are more likely to show large 
angular velocities. Proper motion can be combined with ap- 
parent magnitude to obtain a statistic called the reduced 
proper motion H, which provides a crude estimate of the 
absolute magnitude. 



= m + 5 logio ^1 + 5 

= M + 5 logio Vt - 3.3S 



(2) 



Table 4. Kinematic quantities adopted in this work. The usual 
Galactic coordinate axes UVW are used, with U pointing towards 
the Galactic centre, V in the direction of rotation, and W towards 
the NGP. The velocity dispersion tensor is assumed diagonal in 
this frame. Mean motions are relative to the Sun. 



Although useless for obtaining accurate stellar distances, 
H is sufficient to distinguish populations of stars with dis- 
tinctly different luminosity calibrations or kinematic prop- 
erties. The classical tool for exploiting this is the reduced 
proper motion diagram (RPMD), which plots colour against 
H. The RPMD is topologically equivalent to the HR dia- 
gram, though with considerable vertical scatter due to the 
weak correlation between H and Af . At around ten magni- 
tudes fainter than main sequence stars of the same colour, 
white dwarfs are ideally suited t o identification ba sed on H, 
and several independent studies l|Kilic et al.ll2006l . e.g.) have 
proved this to be a good way to compile clean white dwarf 
samples. 



3.1.1 Tangential velocity selection 

Equation [2] suggests that with an appropriate colour- 
magnitude relation, regions of the RPMD inhabited by white 
dwarfs of different tangential velocity can be isolated. This 
allows us to perform rigorous selections on H to produce 
catalogues of white dwarf candidates within a well defined 
tangential velocity range. This is desirable because cool, low 
velocity white dwarfs overlap in the RPMD with high ve- 
locity subdwarfs from the Galactic halo. Contamination by 
subdwarfs can be reduced by applying a minimum tangen- 
tial velocity threshold to white dwarf candidates, producing 
a cleaner sample of white dwarfs by restricting selection to 
regions of the RPMD more widely separated from the sub- 
dwarf locus. Figure [7] demonstrates the selection of white 
dwarf candidates based on reduced proper motion. The fact 
that low velocity white dwarfs are lost from the survey is 
of course a drawback of this technique; however, the frac- 
tion of stars that fall below the chosen threshold can be 
calculated, if the kinematic properties of the population are 
known. This is done in each field by projecting the veloc- 
ity ellipsoid onto the tangent plane, correcting for the mean 
motion relative to the Sun, and marginalising over the po- 
sition angleto obtain the distribution in tangential velocity 
- see I Murray! ([l9S 3). The values adopted for the mean re- 
fiex motions and velocity disp ersion tensors a r e give n in Ta- 
ble H These are obtained by IChiba fc Beer^ tOOdi) for the 
thick disk and halo; for the thin disk we use the Fuchs et al.l 
(j2009l l study of SDSS M dwarfs, with values taken from 
their 0-lOOpc bin that is least affected by the problems as- 
sociated with the deprojection of pro per motions away from 
the plane (jMcMillan fc Binnevllioogh . Mean motions are rel- 
ative to the Sun; the usual Galactic frame is used in which 
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- ''59F 

Figure 7. Reduced proper motion diagram for all stars in the 
SSA that pass the selection criteria described in Section [2] Over- 
plotted are cooling tracks for H (black) atmosphere white dwarfs 
with vtan = 20 and 40kms~^. Also shown is the cooling track for 
He (grey) atmosphere white dwarfs with vtan = 20kms-i. The 
H tracks are used to identify white dwarf candidates; all objects 
lying below a certain line are selected. 



4 




Figure 8. Two colour plot for all RPM selected WD candidates. 
Squares indicate stars that pass the < 5.0 cut on the model 
atmosphere fit; crosses indicate failures. 



3.2 Photometric parallaxes 

Photometric distances are obtained by fitting the two-colour 
photometry to the white dwa rf model a tmospheres and 
cooling seq uences describe d in Fontaine e t al. (2001) and 
updated in [Bergeron et al.l l|200ll ) and references thereir0. 
They were provided in the SuperCOSMOS bands by 
Dr. Bergeron on request. The models consist of cooling se- 
quences for white dwarfs of different surface gravity and 
H/He atmosphere type. The gravity and atmosphere effect 
the fitted distances by changing the absolute magnitude at 
a given colour, but with only two data points each we can- 
not fit these for our stars. Instead, we assume \ogg — 8.0 
for all our stars, and fit both hydrogen and helium at- 
mospheres. It is well known that the gravities of white 
dwa rfs are tightly distribu ted about this value (for exam- 
ple, [Bergeron etai] (|200ll ) find (logp) = 8.070 ± 0.014), a 
consequence of their narrow mass distribution. Low and high 
mass white dwar fs exist in rou ghly equal numbers (~ 10% 
and 15%) (iLiebert et al.ll2005i ). and fitting to logg = 8.0 
models has opposite effects on the photometric parallax and 
luminosity function. 

Also, the H/He atmosphere type has very little effect 
on the luminosity above around 6, OOOK {bj — tsqf ~ 0.8). 
Below this, the assumption of a H atmosphere for a He at- 
mosphere star will cause the absolute magnitude to be con- 
siderably overestimated, and the distance underestimated. 
In general, optical spectra are useless for distinguishing the 
atmosphere type, because below around SOOO-ff all the ab- 
sorption lines are washed out. There is therefore some am- 
biguity over the nature of the coolest white dwarfs in our 
survey. 



See www . astro . umontreal . ca/~bergeron/CoolingModels 



3. 2. 1 Fitting procedure 

The best fitting logp = 8.0 H and He atmospheres are 
found in a straightforward manner by variance-weighted 
least squares, after interpolating the models at IQK inter- 
vals. The models corresponding to the upper and lower la 
confidence boundaries are found by Xio- = Xmin + 1- We 
take no account of redenning, and do not expect it have a 
significant effect due to the proximity of our stars. Objects 
with Xmin > 5 are rejected from the survey; these are mostly 
unresolved binaries as explained at the end of this section. 

Overall distances are estimated by taking a minimum- 
variance combination of the estimates from each photomet- 
ric band hj, r^gp and ijv. Note that we avoid rn^p/iQ^aS 
in this calculation. Uncertainty on the overall distance es- 
timate is assigned by averaging the upper and lower confi- 
dence boundaries. A two-colour diagram showing the loca- 
tion and status of our stars relative to the models is pre- 
sented in Figure |8l 



3.2.2 Calibration of photometric parallaxes 

In order to estimate the accuracy of our photometric par- 
allaxes, we compare the photometric distances to those ob- 
tained by trigonometri c para llax for a subset of our white 
dwarfs. [Bergeron et al.l (1200 ll ) presented an analysis of 152 
cool white dwarfs with accurate trigonometric parallaxes 
(— < 0.3), all but two of which fall within the proper mo- 
tion and magnitude range of our high proper motion data. 
Of the remaining 150, 116 have SSS counterparts present 
in our input catalogues. We investigated the missing stars 
by searching individual plate records for proper-motion cor- 
rected positions; mostly stars are excluded due to lack of 
a detection at resp/io^aE- Note that when cross identifying 
stars, we do not apply the magnitude and proper motion 
completeness limits, as the identification is done manually 
and contamination is not a problem. We do, however, apply 
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Figure 9. Wc compare the distances obtained by our pliotomct- 
ric parallaxes to those obtained by trigonometric means, for a 
subsample of stars analysed by .Bergeron et al.i (i2001, ). 



the usual restrictions on blend and quality number, as these 
directly affect the quality of the photometry. We also require 
stars to be detected at all four epochs, which is necessary 
for the fit. These constraints reduce the sample to 75 stars. 

We fit atmosphere models to the remaining stars, in 
each cas e using the appropriat e H/He atmosphere as mea- 
sured by [Bergeron et ahl (|200ll ). We place the same cut on 
the residuals as is used in the main survey, in order to re- 
move any stars with spurious photometry. This results in a 
sample of 68 stars with adequately fitted photometric dis- 
tances. A comparison of the distances obtained by the two 
methods is presented in Figure [§1 The correlation between 
the two is r = 0.74, and = 1.08 ± 0.54. The error 

<J-phot 

ir. ^ 



is likely to be lower than this, due to uncertainty 
in dtrig, and we estimate the accuracy of our photometric 
parallaxes crdpi^^t to be around 50%. 

We compared the results of fitting {bj — r^gp ,r5gF — tN) 
and {bj — r^gpjbj — in) colours to the models. The (fej — 
r59F, rsgp — ijv) colours resulted in a slightly worse fit (r = 
0.76, f^"" — 1.09 ± 0.55) presumably due to the superior 

"■photo 

quality of bj, so we adopt {bj — r^gp, bj — iM) for performing 
our photometric parallax fits. We also tried relaxing the 
cut to 6 then 7; in both cases the sample was increased to 
70 stars with = 1.08 ± 0.54. Therefore, aU but a few 

percent of white dwarfs with reliable photometry will pass 
the survey < 5.0 threshold. 



3.2.3 Unresolved binaries 

Selection of white dwarf candidates is made on bj, r^gp and 
/i, so objects with unusual ijv magnitudes turn up at this 
point when their bj In colours are compared to the mod- 
els. The large populations of objects at red bj — in in Figure 
[8] show excess flux in ijv, due to an unresolved cool main 
sequence companion. These WD+dM binaries are ejected 



from the survey by the restriction on the photometric par- 
allax residuals, as with at least one of three bands polluted 
no reliable white dwarf fit is possible. This therefore rep- 
resents a source of incompleteness in the survey. Similarly, 
Sirius-Iike systems consisting of a hot main sequence star 
and white dwarf will fail reduced proper motion selection as 
white dwarf candidates. 

Unresolved double degenerates often have combined 
spectral energy distributions that closely resemble single 
stars of intermediate temperature. They will therefore be 
successfully fitted to the models and pass into the white 
dwarf catalogue as apparently single objects. However, with 
two stars contributing to the fiux the derived photometric 
distance and tangential velocity will be underestimated by 
up to a factor \/2- 

Indications from the local (20pc) white dwarf popula- 
tion are that ^ 70 % of white dwarfs exist as single ob- 
jects iHolberdliopgl ). Around 20% are members of either 
a WD-l-dM or Sirius-like binary system, and so ejection of 
these may result in incompleteness of up to 20%, depending 
on what fraction are unresolved by SuperCOSMOS. Around 
10% of white dwarfs exist in double degenerate binaries, 
which translates to a double degenerate-to-single star frac- 
tion of up to 7% among our catalogue objects, again de- 
pending on what fraction are not spatially resolved. 



3.2.4 Known uHracool white dwarfs 

Of the dozen or so ultracool {T^s < 4000iS') white 
dwarfs that have been reported in the literature, seven 



3250 (iHarris et al. 


I999I; 


(jRuiz & Bereeron 


200 ih. 



l2004h . SSSJ1556 (iRoweU et al, .20081 SDSSJ1452-h45 and 
SDSSJ1632-I-24 (|Harris et al.ll2008h . Our default photomet- 
ric parallax method is inappropriate for these stars for two 
reasons. Firstly, although detailed analysis of stars of this 
class is currently, and necessarily, restricted to only a cou- 
ple of examples, it is clear that their properties are quite 
different to what one would expect based on an extrap- 
olation from higher temperatures. In particular, it would 
appear that most of these objects have extremely He rich 
atmospheres, which is difficult to reconcile with the ex- 
pected accretion rates of H from the interstellar medium 
(jBergeron fc LeggettI [20o3 ) . Also, the single object with a 
trigonometric parallax (LHS 3250) appears over-luminous 
for its temperature, and has been interpeted either as an 
unresolved double degenerate or an extremely low mass sin- 
gle white dwarf. Secondly, models fail to reproduce the SEDs 
of these objects for any set of parameters, indicatin g incom- 
plete input physics. To quote [Harris et all l|2008l ). writing 
with reference to their own objects but applicable more gen- 
erally, 

"It is premature to add these new ultracool white dwarfs to any 
analysis of the space density and luminosity function of white 
dwarfs for two reasons: wc do not yet have models to fit the spec- 
tra adequately to give accurate temperatures and H/He abun- 
dances, and we do not yet have distances to get luminosities, 
masses, and ages." 
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Table 5. Distances and bolometric magnitudes for ultracool 
white dwarfs appearing in our survey. 



Star 


d (pc) 


Vtan (kms ^) 




Tcff{K) 


CE 51 


14.7" 


44 


17.5 


2730 


LHS 3250 


30.3'' 


80 


16.17 


<4000 


LHS 1402 


25° 


58 


16.8 


3240 


SSSJ1556 


32= 


63 


16.17 


<4000 


J1632+24 


23'' 


38 


17.5 


<3000 


SDSSJ0947 


47^ 


18 


16.17 


<4000 


J1452+45 


57'= 


28 


16.17 


<4000 



"Photo tt; 
LHS 3250 ^ 



While this may be true, these stars are present in our survey 
and must be dealt with in some way. With these caveats, we 
proceed to estimate distances and luminosities for these ob- 
jects, though note that our adopted values should be treated 
with caution. 

Of the seven white dwarfs, only three have anything 
close to a reliable distance estimate. LHS 3250 has a trigono- 
metric parallax as noted above, and CE 51 and LHS 1402 
have photometric parallaxes based on spectroscopy and 
multiband photometry, extending into the IR in the latter 
case . SSSJ1556 ha s a SE P very similar to that of LHS 3250, 
and I Rowell et al.l (|2008l ) invoked their similarity to assign 
a distance by assuming these stars also share identical lu- 
minosities. (G^eT^eT^alJIiool) performed a similar analysis 
for SDSSJ0947 using superior Sloan photometry. We con- 
tinue in this way for the final two stars. SDSSJ1452-(-45 is 
closest in colour to LHS3250, though rsg^ — ijv differs by 
~ 0.4 and this star is most likely warmer. SDSSJ1632-I-24 
has identical colours to CE 51 (to ~ 0.01m), and we use 
this star as a reference in this case. The distances and bolo- 
metric magnitudes adopted for these objects are presented 
in Table [5] For all objects, we adopt an identical bolometric 
magnitude error of 0.5M. We have also used SSS proper mo- 
tions to calculate corresponding tangential velocities. Note 
that on the basis of this distance estimate, SDSSJ0947 has 
Vtan = 18kms~^ and drops out of the sample. To be clear, it 
passes the vtan > 20kms~^ RPM threshold for the survey, 
but when the extra (dubious) distance information is added 
at this stage it falls below the cut. 



3.3 Cool white dwarf atmosphere types 

At colours where the choice of atmosphere has a signifi- 
cant effect on the absolute magnitude, a dichotomy arises in 
the photometric d istance estimate that must be addressed. 
iKnox et al.l (|l999l ) deal with the unknown atmosphere type 
by assigning half of the stars H atmospheres and half He. 
However, as He WDs are brighter at a given colour they 
will be sampled over a larger volume and are expected to be 
present in greater numbers than a simple 50:50 ratio. HOG 
use this fact to estimate the relative numbers of the two 
types they expect to find in their survey, in several bolo- 
metric magnitude bins. They also avoid strict atmosphere 
assignments for each object, choosing instead to attach a 
weight to each atmosphere and allow stars to contribute as 
both types. 

We follow the lead of H06, and assign weights to each 



star depending on colour. However, instead of using a few 
discrete magnitude ranges, we derive a continuous weight 
function based on the survey volume for each each type 
as a function of colour. The relative fraction of He to H 
stars at a given bj — r^gp colour is estimated from the cor- 
responding colour-magnitude relations Rnipj — tsqf) and 
RHe{bj — r^gp) assuming a spherical survey volume and 
uniform density profile. We obtain the following formulae 
for the weights uh and uhb for each type; 



riHe 



riHe + riH 



(3) 



where C is the ratio -^^^ of He to H atmosphere white dwarfs 
in a fixed volume. C is likely an evolv ing function of colour, 
due to spectral evolution (see e.g. iTremblav fc BergeronI 
'2008^. Here w e simply set it equal to 0.5 , which is the 
value found bv ITremblav fc BergeronI (|2008l ) for the coolest 
(Tcfi < 10, OOOK) stars in their sample. This covers en- 
tirely the colour range over which the two types differ in 
the colour-magnitude plane, so is an appropriate assump- 
tion for all our stars. The helium weight as a function of 
colour is plotted in Figure [lO(b)| for several assumed values 
of C, along with reference weights obtained if H and He white 
dwarfs did not diverge in absolute magnitude. The diverging 
colour-magnitude relations used to calculate the weights are 
shown in Figure 10(a) 



4 CATALOGUE SUMMARY 

We have undertaken a proper motion and magnitude limited 
survey for white dwarfs in the SuperCOSMOS Sky Survey, 
using two distinct datasets of low and high proper motion 
stars. In both proper motion ranges, most WD candidates 
are ejected by the rio3aB/63F magnitude limit. While this 
is undoutedly the largest restriction on the catalogue size, 
a comparable fraction of low proper motion WD candidates 
are ejected by the low proper motion limits. In light of Sec- 
tion [3]2]2] the catalogue may be up to 50% incomplete due 
to exclusion of blended objects and those for which no first 
epoch detection exists. However, in Section [4.11 we present 
evidence that the incompleteness is uniform within the sur- 
vey volume and does not bias the catalogue. On applying 
a Vtan > 20kms~^ cut in reduced proper motion, we obtain 
9,749 white dwarf candidates with photometric parallaxes 
accurate to around 50% 0. Increasing the Vtan threshold 
results in a cleaner catalogue; vtan > 30kms~^ gives 8206 
stars, and Vtan > 40kms~^ gives 6592. Note that these num- 
bers are based on reduced proper motion selected samples; 
later, we will draw velocity subsamples using the photomet- 
ric parallaxes to determine velocities, and the numbers will 
be slightly different. A sky projection of all WD candidates 
that pass the Vtan > 20kms~^ cut is shown in Figure [TTI 



The catalogue is available for download 
http : //surveys . roe . ac . uk/ssa/links . html 



from 
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Figure 10. (a) Colour-magnitude relations for 0.6Mq H and He 
white dwarf models, showing the divergence beyond bj — r^gp ~ 
0.8. (b) Weights calculated for cool helium atmosphere white 
dwarfs, under different assumptions for the ratio lifis- in the solar 
neighbourhood. The curves correspond to C = 0.1,0.5,2.0 (bot- 
tom to top), and the straight lines are the weights if the H and 
He colour-magnitude relations did not diverge. 



Figure 11. Hammer- Aitoff projection of WD candidates in the 
proper motion survey. Galactic coordinates are plotted with I = 
6 = in the centre. 
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Figure 12. 

20 kms~^ and 30 kms 



distribution for WD candidates with vtan > 
^. The volume calculation assumes the 
density profile and velocity ellipsoid of the thin disk for all stars. 
(y^) = 0.497(0.496) ±0.003 for vtan > 20 (30) kms'i, with 
— following a uniform distribution as expected. Our WD cat- 
alogue is consistent with having been drawn from an unbiased 
sample. 



4.1 Survey bias and incompleteness 

An unbiased survey is uniformly sensitive within the acces- 
sible survey volume. The degree to which our survey meets 
this requirement can be checked by looking at the distribu- 
tion of — , which follows a uniform distribution in the 
ideal case. V is the survey volume contained within the dis- 
tance at which the star lies, and Vmax is the maximum sur- 
vey volume in which the star could reside and still be acces- 
sible to the survey. The calculation of the survey volume is 
non-trivial and will be explained in later Sections; here we 
simply present the distribution, in Figure 1121 The fact the 
distribution is extremely flat is reassuring; this is evidence 
that the large incompleteness present is uniform within the 
survey volume, and does not bias the survey towards any 
particular type of star. As the stellar density profile and 



velocity ellipsoid are included in the calculation of y— — , 
this also suggests that reasonable values have been adopted. 
As -y^ — is drawn from U[0, 1], its distribution should have 
the property that {y^ — ) = 0.5 ± ^^^^ for objects. Our 
Vtan > 20 kms"^ WD catalogue has {y^} = 0.497±0.003. 



4.2 Distance and tangential velocity distributions 

In Figure [TH] we present the tangential velocity distribution 
for all objects in our WD catalogue. Clearly, a small frac- 
tion of stars show velocities much larger than those normally 
associated with the thin disk. Figure [HI shows the distribu- 
tion of distance for all objects in our catalogue. The vast 
majority of stars are within 300pc, suggesting that the lack 
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Figure 13. Tangential velocity distribution of white dwarfs in 
our proper motion survey. 
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Figure 14. Distance distribution of white dwarfs in our proper 
motion survey. 



of a reddening correction does not significantly affect tiie 
photometric parallaxes. 



in the H06 catalogue. Approximately 800 of our WD can- 
didates do not appear, presumably due to slight variations 
in the magnitude completeness limits between USNO-B and 
the SSS re-reductions of the plate material. The WD can- 
didates that are common to both surveys show reason- 
able agreement on stellar parameters: photometric distances 
agree to 85%, and bolometric magnitudes show a disper- 
sion of around 0.5 (assuming hydrogen atmospheres in both 
cases). 



5 LUMINOSITY FUNCTION METHODS 



5.1 The 



density estimator 



There exists a variety of statistical methods for estimating 
luminosity functions, including both parametric and non- 
parametric methods, maximum likelihood estimators and 
simple number counts. The — technique ()Schmidtlll968l ') 
has been used in every major study of the white dwarf lumi- 
nosity function, due largely to tradition and to the relative 
simplicity of the approach, although it has the advantage of 
easily incorporating proper motion selected samples, as well 
as the non-uniform distribution of objects within the Galac- 
tic disk. Its performanc e alongside other methods has been 
analysed bv lGeiio et al.l |2006) and found to be satisfactory; 
it provides an unbiased estimate of the true density, and for 
large enough samples 300) accurately characterises both 
the rising slope and faint peak of the luminosity function. 



The 



method obtains an estimate for the number 



density of objects 4> by summing the inverse of the maxi- 
mum volume in which each object could reside and still be 
accessible to the survey. 



E 



1 



Uncertainties are conventionally assigned assuming Poisson 
statistics, where the standard error in each rr^ — term is 



equal to the term itself (e.g. 1 ± 1 events). These are then 
summed in quadrature to obtain the error on (f), 



2 



E 



1/2 



4.3 Comparison to H06 WD catalogue 

We have compared our WD catalogue to that of H06, which 
was derived from a combination of SDSS DR3 and USNO- 
B data using much the same survey procedure as we have 
applied here. 1927 of our WD candidates with vtan > 30 
kms~^ lie in the DR3 footprint area, compared to 6009 found 
by HOG. The difference is much larger than is accounted 
for by our ~ 50% incompleteness, and is a result of the 
wider survey limits used by HOG. While our surveys reach a 
similar magnitude limit (they are after all based on the same 
photographic plate material), HOG has a much lower proper 
motion limit of /i ^ 0'.'02 yr~^. Our survey is less sensitive 
to proper motion but covers a much wider area, hence the 
comparable catalogue size overall. 

Of the 1927 stars, 1135 have an unambiguous match 



A more accurate approach would be to use the iGehrelsl 
( 1986) upper and lower confidence limits for Poisson statis- 
tics, which would result in a lj!og3 contribution from each 
star for a 68.27% confidence interval. It is to allow compar- 
ison with other studies that we adopt an uncertainty of ±1 
on the number of stars; all previous measurements of the 
white dwarf luminosity function have used this. 

Objects are binned on bolometric magnitude, and the 
density associated with each bin is calcul ated in this manne r 
to obtain the luminosity function. As in lKnox et al.l l|l999f ). 
we plot the luminosity function points at the mean magni- 
tude of the objects in each bin. This is more realistic, and 
shifts the observed luminosity function slightly in regions 
where the number counts change rapidly with magnitude, 
such as at the downturn. We also assign horizontal error 
bars to each point, calculated by averaging the lower and 
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upper bolometric magnitude errors separately, e.g. 



z — 1 up, I 



N 



with aup,i and aiow.i assigned from the la photometric mod- 
els. 



5.2 Calculating Vmax 

Our survey is limited on both apparent magnitude and 
proper motion. The intrinsic stellar properties, namely the 
absolute magnitude and tangential velocity, lead to restric- 
tions on the distance at which each star could reside and 
still pass the survey limits. The apparent magnitude limits 
are fixed in each field, and the corresponding distance limits 
are found according to 



fO" 



(4) 



f^mm = max f 



where the index i iterates over each of the four bands. 

Due to the non-analytic lower proper motion limits es- 
tablished for each field, no simple expression exists for the 
corresponding survey distance limits. Indeed, if the lower 
limit changes rapidly with apparent magnitude there may 
even be several ranges of distance in which the star passes 
the survey limits. Vmax must be calculated by integrating the 
appropriate stellar density profile ^ along the line of sight 
between cij^i„ and d^ax , at each step evaluating whether the 
star passes the proper motion limits, which are calculated 
from the star's apparent bj magnitude at that distance. This 
leads to the integral 



Vmax ~ y 



V{r)r .dr 



where the summation is over all surveys fields, and 



V{r) = 



1 if 



/^min(^j(^)) 4 74 J. ^ f^max 



otherwise. 



This method for Vmax follows that of IStobie et al. J 19891). 
genera lised to arbitrary Galactic latitudes bv iTinnev et all 
l|l993l ) and further here to allow for the piecewise lower 
proper motion limits. 



5.2.1 Stellar density profiles 

For the thin and thick disks we use an exponential decay law 
in Galactic plane distance |z»| for the stellar density profile, 
in order to correct for the truncation of the survey volume 
by the scaleheight effect. We expect to see no effect arising 
from the radial scalelength of the disk, due to the relative 
proximity of our stars. The appropriate form for ^ is thus 



— = exp ■ 

P© 



H 



where H is the scaleheight. We adopt H = 250pc for the thin 
disk, w hich is in line with the result of iMendez fc GuzmanI 
l|l998l ) obtained for faint main sequence stars. These are 



likely of similar age to the white dwarfs in our catalogue 
and are expected to show a similar spatial distribution, hav- 
ing been subjected to the same kinematic heating. This is 
also the value used in most other studies of the white dwarf 
luminosity function, and thus allows more meaningful com- 
parison with other works. There is some empirical evidence 
that the scaleheight of disk white dwarfs increases towards 
fainter magnitudes (see HOG) where the stars are on average 
older, but we ignore this here. At all but the brightest mag- 
nitudes our white dwarfs are so close to the Sun that the 
chosen scaleheight makes very little difference anyway. 

The Solar distance from the Gala ctic plane, zq, is o f- 
ten omitted in studies like this (e.g. iTinnev et al.|[l993l ). 
which is equivalent to setting its value to zero. However, the 
consensus of many star count inves tigations is that in fact 
zq lies close to ~ 20pc (lReed|[2"006l ). The adjusted density 
profile becomes: 

P -|r sin(b) -|- zq\ 



= exp ■ 



H 



where h, r are the Galactic latitude and line of sight distance 
and Zq = 20pc is the Galactic plane distance of the Sun. 

For the spheroid, we use a uniform density profile. We 
expect to see no variation in stellar density over the distances 
probed by our survey. 

5.2.2 Corrections 

Several steps in the compilation of our white dwarf catalogue 
have the side effect of excluding a fraction of target stars. 
We correct our density estimate for the ejected objects, un- 
der the assumption that the incompleteness is uniform with 
luminosity and does not bias the survey towards any partic- 
ular type of star. The discovery fraction x of stars that pass 
the tangential velocity threshold is calculated as described 
in Section 13.1.11 and the contribution of each star to the 
total density is adjusted according to 



^ 1 



-j' X* Vmax.i 

where x for each star is taken from the field in which the 
star was discovered. A similar correction arises from the re- 
striction on astrometric residuals described in Section 12.1.31 

5.2.3 Atmosphere types 

Fainter than Mboi = 14 {bj — rsgp ~ 0.8), the choice of H or 
He atmosphere solution has a significant effect on the fitted 
bolometric magnitude. In order to account for the unknown 
H/He atmosphere types of cool white dwarfs, stars are al- 
lowed to contribute as both types, with a weight set accord- 
ing to Equation [3l In each case, — is simply multiplied 
by the appropriate weight before inclusion in the sum: 



H 

max A 



He 

max A 



Note that Vmax is different for each solution, and they will 
not in general contribute to the same luminosity bin. Also, 
the different photometric distances lead to different tangen- 
tial velocities, and in some cases only one of the two solutions 
will pass the velocity threshold and be included in the LF. 
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THE WHITE DWARF LUMINOSITY 
FUNCTION 



In Figure 15(a) we present the luminosity function for white 
dwarfs in the SSS, on adopting a 250pc scaleheight and min- 
imum tangential velocity threshold of 30kms~^. The struc- 
ture in the luminosity function at the faint end is easily 
discernible - beyond the peak, there is a sharp drop off fol- 
lowed by a slow decline. Theory predicts that high mass 
white dwarfs cool faster than their normal mass counter- 
parts, and, all other things being equal, fall in the region 
beyond the peak where the luminosity function for normal 
mass white dwarfs terminates. Good constraint in this re- 
gion is vital for obtaining accurate age estimates, and the 
number of datapoints beyond the peak is encouraging. How- 
ever, a quantitative analysis is only possible in conjunction 
with theoretical luminosity functions. 

Due to the effect of the magnitude-dependent proper 
motion limits, at fainter bolometric magnitudes the sam- 
ple is dominated by stars of brighter apparent magnitude. 
This improves the photometric parallax fit and reduces the 
width of the horizontal error bars on the luminosity func- 
tion points. Beyond the peak, the sample is dominated by 
the ultracool white dwarfs, which have a fixed bolometric 
magnitude uncertainty of 0.5M leading to larger error bars 
relative to the peak. 

Figure |15(b)| shows the effect on the luminosity func- 
tion when the lower tangential velocity threshold is varied. 
The fact that no systematic trends are apparent in the lumi- 
nosity function suggests that the sample is not significantly 
contaminated by subdwarfs at low velocities. 



6.1 The luminosity function for high velocity 
white dwarfs 

Any spheroid white dwarfs present in our catalogue may 
be identified by their large tangential velocities. Figure [T6l 
shows the tangential velocity distributions for the thin disk, 
thick disk and spheroid along the line of sight to one of our 
survey fields, as determined from M dwarfs and low m etallic- 
ity stars fsee lFuchs et al.ll2009l : [chiba fc Beersll2000l ). A cut 
of Mt > 200kms~^ is often considered to cleanly separate the 
spheroid and disk populations, and the luminosity function 
obtained on applying this cut to our catalogue is presented 



in Figure 17(a) In this case, the discovery fractions used 



to correct for the excluded low velocity stars are calculated 
from the spheroid velocity ellipsoid, and the density pro- 
file is that of a uniform population. The effect of varying 
the velocity threshold is investigated in Figure 17(b) The 
fact that the vt > 160kms~^ luminosity function sits at a 
slightly higher density suggests that there is some residual 
contamination from the disk even at these velocities. 

The origin of the high velocity white dwarf component 
can be further probed by looking at the bulk motion. At- 
tributing individual stars to a particular kinematic popu- 
lation based on tangential velocity alone is tricky, because 
the kinematic properties of cool white dwarfs are relatively 
uncertain, and it is possible that the high velocity tail of the 
disk population (s) overlaps considerably with the spheroid 
(see iReidI l2005f ). Ideally, radial velocities would complete 
the full 3D space motion in Galactic coordinates, allow- 
ing far better discrimination for individual stars. However, 
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Figure 15. (a) Luminosity function for vt > 30kms~^ white 
dwarfs in the SSS. (b) Luminosity functions for white dwarfs on 
adopting a range of lower tangential velocity cuts. 
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Figure 16. Tangential velocity distributions for the three major 
kinematic populations, along the line of sight to field 362 in the 
southern hemisphere. A cut of vt > 200kms~^ is often considered 
to cleanly separate the spheroid and disk populations. 
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Figure 17. Luminosity functions for high tangential velocity 
white dwarfs. Figure (a) shows the LF for vt > 200kms~^ WDs, 
and figure (b) compares the LFs obtained when the lower velocity 
limit is varied. 



with only proper motions it is still possible to measure the 
mean motion in Galactic coordinates for the population as 
a whole, which is sufficient to distinguish a spheroid sam- 
ple from one drawn from a rotating disk. This is done by 
deproje cting the proper motions , according to the method 
used bv lDehnen fc BinnevI (|l998l ) to analyse the kinematics 
of stars in the Hipparcos catalogue. 

On doing so, we find that the 93 stars with Vt > 
200kms-i have (UVW) = (-88 ± 12, -214 ± 24, -38 ± 5) 
kms"'^. These would suggest that the high velocity members 
of our WD catalogue samples are drawn from a non-rotating 
population, i.e. that of the spheroid. We therefore conclude 

is 



that the luminosity function presented in Figure 17(a 
representative of the spheroid white dwarf population. 



6.2 Local disk and spheroid white dwarf densities 

Integrating the luminosity function presented in Figure 



15(a) gives a total local density for white dwarfs in the solar 
neighbourhood of (3.19 ± 0.09) x 10"^pc"^ 

The spheroid white dwarf luminosity function of Fig- 
T7(a)] integrates out to (4.4 ± 1.3) x 10~^pc~^, resulting 



in a disk-to-spheroid ratio of 725 ± 215 disk stars to every 
spheroid star. These densities are significantly lower than 



thos e obtained by o ther s tudies, e.g. H06 find 4.6 x 10 pc 
and iLeggett et all (|l998l ') find 3.4 x 10"^pc"^ for the local 
density of disk white dwarfs, and H06 find 4 x 10^^ pc^"^ for 
those of the spheroid. If we recall from Sections l3.2.2l and l4.1l 
our catalogues may be up to 50% incomplete due to blended 
objects and those missed at first epoch r, which explains the 
disagreement between these numbers. This incompleteness 
is expected to be uniform with bolometric magnitude, and 
therefore will not affect the disk-to-spheroid ratio, nor will 
it affect any conclusions about the age of these two pop- 
ulations, as this is insensitive to the normalisation of the 
luminosity function. It should be noted at this point that 
while the disk WDLF almost certainly falls rapidly to zero 
beyond the magnitude range of our survey, the spheroid 
WDLF most likely continues rising to fainter magnitudes. 
Deeper surveys may find significantly larger integrated den- 
sities for the spheroid by including fainter magnitude bins, 
and the numbers presented here should be regarded as lower 
limits. The contribution from fainter objects may also be 
constrained using theoretical luminosity functions, however 
we do not do this here. 



7 UNTANGLING THE DISKS AND SPHEROID 

The technique of drawing velocity selected sub-samples of 
stars works reasonably well at isolating a clean sample of 
spheroid stars. Although uncontaminated by disk stars, sam- 
ples drawn in this manner do of course miss a large fraction 
of spheroid stars at lower velocities, where they are indis- 
tinguishable from the disk. Also, the thin and thick disks 
could never be separated in this manner, due to the large 
overlap in tangential velocity leaving no range in which the 
populations can be reliably isolated from each other or the 
spheroid. It is worth pointing out that chemical tagging, as 
might be used to distinguish the main sequence members of 
these populations, cannot be used for white dwarfs as any 
photospheric metals sink rapidly below the photosphere. 

An alternative approach can be derived by considering 
the how the sampled volume for each of the kinematic pop- 
ulations affects the numbers of stars that make their way 
into the survey. The total number of survey stars, A^, , is 
determined by the local number density n multiplied by the 
survey volume V , separately for each kinematic population, 

N, = nthin X Khin + "-thick X Vthick + "-sph. X Vsph. (5) 

where the subscripts refer to each of the major kinematic 
populations. In a narrow range of magnitude, such as one 
of our luminosity function bins, the factor that determines 
V for each of the populations is the tangential velocity dis- 
tribution combined with the survey proper motion and tan- 
gential velocity limits. The magnitude limits are effectively 
decoupled from the analysis, because stars belonging to each 
population have very similar mean absolute magnitudes over 
the small ranges considered. The kinematics differ consider- 
ably however, and for equal survey limits each population 
will be sampled over a different volume of space. By vary- 
ing the tangential velocity limit, and recalculating each V 
for the new Nt, a set of linear equations in the unknowns 
densities n can be generated. In general, the equation set is 
non-singular and solvable using linear algebra techniques. 
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The power of this approach Hes in the fact that individ- 
ual stars are not assigned conclusively to either population. 
Instead, it simply measures the fraction of stars that belong 
to each population, as a function of vt- Note that we will re- 
fer to the survey volume measured in this technique as the 
effective volume, Veff, to distinguish it from that used in the 
conventional y-^ — technique. In this approach no strict dis- 
tance boundaries arise from the proper motion limits and 
the survey is sensitive to some velocity sub-sample at all 
distances. 



second argument in each case is the velocity limit arising 
from the survey proper motion limits at distance r. 

The velocity ellipsoids adopted are those of Table [l] 
As before, a uniform density profile is used for the spheroid 
and a 250pc scaleheight exponential disk is used for the thin 
disk. The thick disk density profile is that of a 1500pc scale- 
height exponential disk; however, the proximity of our stars 
means the effective volume is rather insensitive to the scale- 
height at this level, and it is the velocity correction alone 
that separates the thick disk and spheroid populations in 
our survey. 



7.1 Modelling the effective volume 

The survey volumes used in this analysis are calculated in 
a rather different manner to Vmax from Section [5.21 This is 
because neither the absolute magnitude nor tangential ve- 
locity are directly observed from the survey objects. In the 
first instance, we instead use the theoretical mean absolute 
magnitude for all stars in each luminosity bin, obtained by 
integrating model M{Mboi) relations for each survey band. 
We assume the LF is flat over the width of each bin, and 
use standard logg = 8.0 DA models; the difference between 
the effective volumes for DA and DB models is very small 
at fixed bolometric magnitude. These are then used in con- 
junction with the apparent magnitude limits in each survey 
field to place limits on the distance at which this hypothet- 
ical star could lie and still pass the magnitude limits. These 
are found according to 

'"i.max-<Mi> 



max I 10 



in 



which is similar to the expression in Equation [4] except the 
observed absolute magnitude has been replaced by the the- 
oretical mean for the LF bin. 

The effective survey volume for each kinematic popu- 
lation is then found by integrating the appropriate density 
profile between d™i„ and d™aa;, at each step correcting 
for the fraction of stars that pass the tangential velocity lim- 
its at that distance. The survey volume is thus marginalised 
over the tangential velocity, leading to the effective volume 
probed for each population. This is done as follows. 



where the summation is over all survey fields. The factor 
is the discovery fraction of stars that pass the tangen- 
tial velocity limits at distance r, and is calculated from the 
cumulative distribution by 

X{r) = cdi{{Vupper{r)) - cdf{{vionier{r)) 

and viower are tangential velocity limits fixed 



where Vu 

by the survey design. The appropriate Vupper 
found according to 



and Viewer are 



Vupper (r) = rain{vmax , 4.74 fj,max r) 
viower{r) = max(iimi„,4.74^min r) 

where Vmax and Vmin are external tangential velocity limits 
that may be applied to restrict the velocity range, and the 



7.2 Solution for the number density 

Applying external Vmin and Vmax cuts enables one to gen- 
erate several instances of Equation O which can be used to 
construct a linear equation set in the unknown densities. We 
use non-overlapping vtan ranges to construct each equation, 
so that the covariance of neighbouring star count bins is 
zero. The set of equations can be cast in matrix form like 



Vthin{5vi) Vthick{5vi) Vsph.(5vi) 
Vthiu{5v2) Vthick{Sv2) VsphXSv2) 
Vthin{5V'j,) Vthick{5vz) VsphX^Vz) 



VthiniSVm) Vthick{SVm.) VsphX^Vm). 



N,{Sv2) 

where Svi^2 m are the chosen tangential velocity ranges. 

Investigations indicate that a suitable set of velocity ranges 
is 30 < -ui < 50kms-\ 50 < W2 < 80kms~\ 80 < «3 < 
120kms"\ 120 < t;4 < 200kms"^ and 200 < us < ookms"\ 
thus sampling at finer resolution at highly populated ve- 
locities where the relative contributions of each population 
change rapidly, and probing the region of pure spheroid stars 
at the extreme. In short hand, 

Vn^N 

We assume no errors on the design matrix V , permitting a 
reasonable solution to be obtained by weighted least squares. 
The weighted least squares solution for n, denoted fi, is given 
by 

n = {V''WV)-\V''WN) 

where ^ denotes the matrix transpose. We use a non-negative 
least squares algori thm to enforce the posit ivity constraint 
on the parameters (|Lawson &: Hansonlll974f ). W is the ma- 
trix of weights, which we set according to 
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i.e. inverse variance weights, adopting Poisson statistics to 
estimate the noise on the observed number counts. For dis- 
crete tangential velocity ranges, there is no covariance be- 
tween adjacent number count bins and the weight matrix 
is diagonal. The uncertainties on n are obtained from the 
variance-covariance matrix V calculated 

V = (yi"M/v)-' 

and we adopt these as the formal uncertainties on the LF 
points. Note that the use of theroretical mean absolute mag- 
nitudes removes the horizontal error bars. 



7.3 The luminosity functions for the disks and 
spheroid 

The luminosity functions for the thin disk, thick disk and 
spheroid obtained on applying this technique to our white 
dwarf catalogue are presented in Figure 1181 Note that we 
do not connect up luminosity function points across empty 
magnitude bins, where no solution was found for one or more 
of the populations. The turnover in the thin disk is apparent 
at around Mi,ai ~ 15.75, in agreement with other works. In 
contrast, the spheroid luminosity function continues to rise 
to the faintest detected magnitudes. No solution is found 
for the spheroid in a single bright magnitude bin at M(,oj = 
12.75, which is most likely due to statistical fluctuations 
in the spheroid number counts rather than a failure of the 
model. 

The thick disk emerges nicely as an intermediate density 
population between the thin and thick disks. There is some 
indication of a turnover at faint magnitudes, but again the 
constraint is too poor to draw firm conclusions and beyond 
Mtoi = 15.25 this population is lost altogether. The fact the 
thin and thick disks reach approximately equal density at 
magnitudes brighter than Mboi ~ 7 is interesting, and may 
hold information about the relative star formation histories 
of these populations. However, the xt statistic presented in 
the lower right tile indicates that a two component model 
excluding the thick disk would fit the data only marginally 
worse than the three component model at these magnitudes, 
so the detection of the thick disk as a distinct entity is not 
conclusive. It may be the case that the scaleheight or velocity 
ellipsoid of the thick disk white dwarfs is closer to that of 
the thin disk at these magnitudes. 

We applied this technique to our white dwarf catalogue 
several times using different combinations of the kinematic 
populations, to check that the model required all three pop- 
ulations to fit the data. The presence of the spheroid stars 
can be proved by conventional techniques (c.f. Section [6. ip . 
but the thick disk can only be detected by looking for a sig- 
nificant improvement in the fit when it is included in the 
model. The upper xt plot indicates that models with just 
a thin disk, or with a thin disk and spheroid, fit the data 
very poorly. In the lower panel we have included the thick 
disk; the fit is clearly much better, so we regard this as the 
first detection of the thick disk white dwarf luminosity func- 
tion. The fit is improved marginally by the inclusion of the 
spheroid population. 

The faintest few bins are very poorly constrained, and 
are significantly more uncertain than the constraint achieved 
using the traditional — method. This is disappointing be- 
cause good constraint at these magnitudes is necessary for 
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Figure 19. Comparison between LFs obtained by traditional 
— and our effective volume technique. In the case of the 
thin disk l |19(a)[ l, the lower density found by the effective vol- 
ume technique is due to the fact thick disk objects have been 
excluded. Both techniques agree well on the spheroid luminosity 
function because a clean sample can be isolated for the — 
method. However, the effective volume technique has a smaller 
error because the spheroid stars at low velocity are included in 
the solution. 



accurate age measurement. Most of the age information is 
contained in the depth and position of the turnover at faint 
magnitudes, although quantitative conclusions require anal- 
ysis in conjunction with theoretical luminosity functions. 
The three component model also fails at these magnitudes as 
no positive solution is found for the thick disk. This means 
the densities recovered for the thin disk and spheroid are 
somewhat compromised. The reason for this may be due to 
the fact the faint bins are dominated by the ultracool white 
dwarfs, which have particularly uncertain photometric par- 
allaxes and bolometric magnitudes. It may also be the case 
that the both the velocity ellipsoid and scaleheight of the 
faintest disk white dwarfs are infiated relative to brighter 
stars, due to the faint stars being on average older. This 
would result in incorrect effective volumes and an inappro- 
priate model. In Figure 1191 we compare the thin disk and 
spheroid luminosity function obtained using the effective 
volume approach to those obtained from the same star cata- 
logue using the conventional — technique. In the case of 
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Thin + thick disks 
Thin disk, thick disk + spheroid 




Figure 18. The luminosity functions for thin disk, thick disk and spheroid white dwarfs, obtained from our catalogue using 
the effective volume method. The lower right tile shows the statistic for each luminosity bin, using different models to 
predict the observed star counts. 



the thin disk ( 19(a) I, the lower density found by the efTec- 



tive volume technique is due to the fact thick disk objects 



have been excluded. The traditional 



technique used in 



previous studies of the disk white dwarf luminosity function 
does not distinguish thin and thick disk stars, so in this case 
and all others the 'thin disk' luminosity function is really 
a sum of both the thin and thick disks. Both techniques 
agree well on the spheroid luminosity function, because a 
clean sample can be isolated for the -ry^ — method. How- 
ever, the effective volume technique has smaller error bars 
on all points because the spheroid stars at low velocity are 
included in the solution. It is also able to probe 2.5 magni- 
tudes deeper, although only one of these bins has reasonable 
constraint on the density. 



7.4 Thin disk/thick disk/spheroid normalisation 

Integrating the luminosity functions in Figure [18] allows us 
to measure the total and relative densities of the different 
kinematic populations in the Solar neighbourhood. For the 
thin disk, we find a total density of (2.23±0.17) x lO'^pc"^. 
This is somewhat smaller than that measured in Section [6^ 



due to the exclusion of thick disk objects her^. For the thick 
disk, we find a total density of (4.55 ± 0.41) x 10~''pc~^, 
and for the spheroid we find (1.4 ± 5.6) x 10"*pc"^. The 
spheroid density is considerably larger than that of Section 
I6.2l due to the fact the luminosity function extends to much 
fainter magnitudes where the density is greater. However 
the same caveat applies; the spheroid WDLF may extend 
to considerably fainter magnitudes than probed here, so the 
density should be regarded as a lower limit. 

We can correct these values for the incompleteness in 
our survey if we assume the difference between the disk den- 
sity measured in Section 16.21 and that measured by HOG 
is due to the incompleteness. This boosts the densities to 



3.1 X 10-^pc-^ 6.4 X lO-^pc" 



^ and 1.9 x lO'^pc"^ for the 



thin disk/thick disk/spheroid Using these numbers, the total 
local density of white dwarfs is split among the kinematic 
populations in roughly 79%/16%/5%. This agrees well with 
predict ions based o n Galactic models, in particular the re- 
sults of lReidI l|2005h who adopts numbers similar to these to 
successfully reproduce the results of a number of proper mo- 
tion surveys for white dwarfs. The similarity between Reid's 



* Note that the thick disk density is not equal to the difference, 
because in Section l6.2l the Vmax for each thick disk object is cal- 
culated using the thin disk density profile and velocity ellipsoid. 
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Figure 9 and our luminosity functions in Figure [18] is strik- 
ing. 



8 COMPARISON TO OTHER WORKS 

The most directly comparable white dwarf luminosity func- 
tion to this study, in terms of number of stars and survey 
technique, is that of H06. These authors used proper mo- 
tions derived from a combination of SDSS and USNO-B 
astrometry to obtain a sample of 6000 white dwarfs with 
Vt > 30kms~^, with photometric parallaxes obtained from 
superior 5 band SDSS photometry. Their disk luminosity 
function is shown in Figure [20 (a) [ alongside our own repro- 
duced fro m Figure | 18[ At t he br ight end, the luminosity 
function of iKrzesinski et al] (|2009l ) is shown, also obtained 
from the SDSS but using selection based on UV excess to 
assemble a sample of hot white dwarfs. The most obvious 
difference in the thin disk luminosity functions (top panel) 
is the rather large vertical offset, which is due to a com- 
bination of the incompleteness in our survey, and the fact 
that our luminosity function is exclusively for thin disk ob- 
jects whereas those shown for comparison contain thick disk 
objects as well. At the bright end the luminosity functions 
diverge, with ours finding a considerably lower density for 
hot white dwarfs. Our luminosity function obtained using 
the T 



method agrees much better (e.g. Figure 19(a) 



so the difference may be due to a genuine change in the 
relative contributions of the thin and thick disks at these lu- 
minosities. However, Section [7.31 suggests the thin and thick 
disk decomposition may be unreliable in this range. There 
are some striking similarities between these luminosity func- 
tions too. In particular, they agree very well on the position 
of the turn over at the faint end, a feature arising from the 
finite age of the Galaxy. Our — luminosity function finds 
the same result, while achieving greater constraint. The in- 
flexion at Mtoi ~ 11 identified by H06 seems to be confirmed 
by this study. This feature may hold valuable information 
about the recent star formation history of the disk. 

The spheroid luminosity functions (lower panel) show 
the same vertical offset, though this is less significant due 
to the larger errors. Morphological features of note include 
the slope of the luminosity function at brighter magnitudes 
where the constraint is good, and the rather more uncertain 
upturn at around Mf,oi ~ 15. Both surveys seem to agree on 
these features. Our spheroid luminosity function has consid- 
erably smaller error bars at magnitudes where the surveys 
overlap, and probes two and a half magnitudes fainter. 



9 CONCLUSIONS 

9.1 Proper motion survey and the luminosity 
function for white dwarfs 

We have conducted a wide-angle survey for white dwarfs us- 
ing data from the SuperCOSMOS Sky Survey and reduced 
proper motion selection. Our catalogue contains 9749 WD 
candidates on applying a Vtan > 20kms~^ threshold, with 
photometric distances accurate to around 50%. Our cata- 
logue suffers from ~ 50% incompleteness due to the ejection 
of blended objects and those undetected at the first epoch, 
however this has been shown to be uniform and does not 
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Figure 20. Figure (a) compares the decomposed thin disk lu- 
minosity function deri ved in this work with that of H06 and 
IKrzesinski et al. who use traditional methods. The ver- 

tical offset is a combination of uniform incompleteness and the 
ejection of thick disk objects from our luminosity function. Fig- 
ure (b) shows the H06 luminosity function for white dwarfs with 
Vt > 200kms~^, as well as our decomposed spheroid luminosity 
function. 



bias the survey. Although the normalisation of the result- 
ing luminosity functions is compromised their morphology 
is not, and any cosmochronology measurements using lumi- 
nosity function models will likewise be unaffected. 

We have introduced a new technique to decompose a 
single mixed catalogue of stars into the contributions from 
each kinematic component, and used this to measure the 
thin disk, thick disk and spheroid WDLFs separately. This 
is the first direct measurement of the thick disk WDLF. Our 
thin disk and spheroid WDLFs probe 1.0 and 2.5 magnitudes 
deeper than the next deepest study (H06). By integrating 
the WDLFs we have measured the relative contributions 
to the total WD density in the Solar neighbourhood from 
each of the kinematic populations, finding results in good 
agreement with predictions. 

The excellent agreement between our thin disk WDLF 
and that of HOB confirms the location of the faint turnover 
at Mtoi = 15.75. This is surely now a secure result. An 
interesting feature in the rising slope of the thin disk WDLF 
that is also observed by HOB would seem to be confirmed by 
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our results. At M^oi ~ 11 WD cooling times are of the order 
500 Myr, so this feature may arise from recent star formation 
activity in the disk. 

9.2 Decomposing the kinematic populations in 
the Solar neighbourhood 

Our work presents a useful way to decompose the otherwise 
indistinguishable WD components of the Galactic kinematic 
populations in the Solar neighbourhood, although our sur- 
vey does not probe faint enough magnitudes to reach the 
WDLF terminus in each case. This technique may be applied 
to other types of star, although it is of particular applica- 
tion in WD studies due to the fact that metallicity cannot 
be used to assign population membership, and the lack of 
spectral lines in cool stars means radial velocities are not 
available. 

As our method separates the populations based on their 
tangential velocity distributions, the fact that there has been 
no rigorous study of the kinematics of faint white dwarfs 
means our velocity dispersions (adopted from main sequence 
stars) must be considered a source of error. For example, as- 
suming a thin disk velocity ellipsoid that is narrower than 
the true value may result in the thick disk being over popu- 
lated in our models. It may also be the case that WDs at the 
faint end of the thin disk luminosity function have undergone 
significantly more kinematic heating than those at brighter 
magnitudes, and no single velocity dispersion is truly rep- 
resentative. Errors on the photometric parallaxes will also 
effect the partitioning of stars among the kinematic popula- 
tions, as they propagate directly into errors on the tangential 
velocity. The magnitude of these effects could be estimated 
using suitable Monte Carlo models, however we have not 
done this here. 

We also point out that while the moments of the veloc- 
ity di stribution can be recover ed from proper motions alone 
(e.g. iDehnen fc Binnevlll998[ ). reduced proper motion se- 
lected catalogues such as our own are most likely useless for 
measuring velocity distributions, due to the fact that the 
catalogue excludes low velocity stars. 

9.3 A note on WDLF age estimates 

The technique of dating a population of stars by measuring 
the position of the WDLF turnover has been applied suc- 
cessfully to young Galactic clusters. However, the position of 
the turnover is relatively insensitive to age at the canonical 
disk age of 8 Gyr, and it is the depth and structure beyond 
the turnover that is the main source of age constraint in 
this regime. This structure is still extremely uncertain for 
several reasons, including statistical error arising from the 
small numbers of stars and uncertain bolometric magnitudes 
for the ultracool white dwarfs that populate this range. Our 
work has shown that the different Galactic kinematic popu- 
lations overlap significantly beyond the thin disk LF peak, 
and traditional approaches to the WDLF that do not allow 
for the possibility of a significant contribution from thick 
disk and spheroid stars cannot hope to measure an accurate 
thin disk age. For these reasons it is our opinion that ex- 
isting Galactic disk age measurements based on the WDLF 
turnover should be treated with caution. 



Accurate aging of the kinematic populations separately 
would indeed be a major achievement of WD cosmochronol- 
ogy. This may have to wait for results from the PanSTARRS 
or Gaia surveys before the observational resources are in 
place. These would have to be matched by an improved un- 
derstanding of the evolution and cooling times of the faintest 
white dwarfs, so that some of the current systematic errors 
(e.g. the spectral evolution problem) can be reduced. 
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